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Abstract 

We show that the Politzer theorem on the equations of motion implies approximate 
constraints on the quark correlator, restricting considerably, for sufficiently large Q^, 
the number of independent distribution functions that characterize the internal 
structure of the nucleon, and of independent fragmentation functions. This result 
leads us to suggesting an alternative method for determining transversity. Moreover 
our approach implies predictions on the Q^- dependence of some azimuthal 
asymmetries, like Sivers, Qiu-Sterman and Collins asymmetry. Lastly, we discuss 
some implications on the Burkhardt-Cottingham and Efremov-Leader-Teryaev sum 

rules. 

PACS numbers: 13.85.Ni, 13.88.+e, 11.15.-q 



1 Introduction 



The problem of calculating inclusive cross sections at high energies and high momen- 
tum transfers has become quite important in the last two decades, during which a lot 
of experimental data on deep inelastic processes have been accumulated. In particular 
we refer to deep inelastic scattering (D1S)[1-12], semi-inclusive DIS (SIDIS) [13-26], 
Drell-Yan (DY) [27-32], e+e~ annihilation into two back-to-back jets [33], while analo- 
gous experiments have been planned recently [34- 39]. One of the aims of high energy 
physicists is to extract from data distribution and/or fragmentation functions, espe- 
cially if unknown. Among them, the transversity [40-43] is of particular interest, since 
it is the only twists-2 distribution function for which very poor information [44, 45] is 
available till now. But also transverse momentum dependent (TMD) functions - es- 
pecially the T-odd ones - are taken in great consideration; for instance, knowledge of 
the Collins fragmentation function [46] or of the Boer-Mulders function [47] could help 
extracting transversity, which is chiral-odd and therefore couples only with chiral- 
odd functions. Moreover, TMD functions are involved in several intriguing azimuthal 
asymmetries, hke the already mentioned Collins[46] and Boer- Mulders [47] effects, or 
the Sivers[48, 49], Qiu-Sterman[50, 51, 52] and Cahn[53, 54] effects, which, in part, 
have found experimental confirmation [19, 20, 21, 22, 23, 25, 33] and, in any case, 
have stimulated a great deal of articles [55-64]. Lastly some questions remain open, 
among which the parton interpretation of the polarized structure function 5'2[65, 66]. 
Obviously, all of these data and kinds of problems are confronted with the QCD the- 
ory and in this comparison short and long distance scales are interested, so that the 
factorization theorems[67-70] play a quite important role in separating the two kinds 
of effects. Strong contributions in this sense have been given by Pohtzer[71], EUis, 
Furmansky and Petronzio[72, 73](EFP), Efremov, Radyushkin and Teryaev[74, 75], 
CoUins, Soper and Sterman[76, 77, 69], and Levelt and Mulders[78](LM). 

In the present paper we propose an approach somewhat similar to EFP's and 
to LM's, but we use more extensively the Politzer theorem on equations of motion 
(E0M[71]). We consider in particular the hadronic tensor for SlDlS, DY and e^e~ — > 
ttttX. We also consider energies and momentum transfers high enough for assuming 
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one photon approximation, but not so large that weak interactions be comparable 
with electromagnetic ones. Lastly, as regards time-like photons, we assume to be far 
from masses of vector resonances, hke J/^', T or 

Our starting point is the "Born" approximation[78] for the hadronic tensor, which 
reads, in the three above mentioned reactions, as 

W^p{Pa. Pb, q)^C^el[ dp^d^p^Tr [niPha^Whp] ■ (1) 

Here C is due to color degree of freedom, C = 1 for SIDIS and 1/3 for DY and e^e" 
annihilation, p and p' denote the four-momenta of the active partons, such that 

= (2) 

q being the four-momentum of the virtual photon and the — sign referring to SIDIS, 
the + to DY or to e^e^ annihilation. $^ and <I>b are correlators, relating the active 
partons to the (initial or final) hadrons Ha and /i^, whose four- momenta arc, respec- 
tively. Pa and Pb- We restrict ourselves to spinless and spin-1/2 hadrons. a and 
h are the flavors of the active partons, with a — u,d, s, u, d, s and 6 = a in SIDIS, 
6 = a in DY and e'^e~ annihilation; is the fractional charge of flavor a. In DY 
and ^B encode information on the active quark and antiquark distributions inside 
the initial hadrons. In SIDIS $b is replaced by the fragmentation correlator A^, 
describing the fragmentation of the struck quark into the final hadron Kb- In the case 
of e+e^ annihilation, both correlators $^ and have to be replaced by and A^ 
respectively. 

In the approximation considered we define the distribution correlator (commonly 
named correlator) as 

$,,(p;P,5) =nJ -0^^e^^^{P,S\^^{Q)A{x)\P,S). (3) 

Here N is a normalization constant, to be determined in sect. 4. ip is the quark* 
field of a given fiavor and \P,S) a state of a hadron (of spin or 1/2) with a given 
four-momentum P and Pauli-Lubanski (PL) four-vector S, while p is the quark four- 
momentum. The color and flavor indices have been omitted in ijj for the sake of 



*For an antiquark eqs. (3) and (4) should be slightly modified, as we shall see in sects. 2 and 6. 
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simplicity and from now on will be forgotten, unless differently stated. On the other 
hand, the fragmentation correlator is defined as 

A,,.(p; P,S)^nJ ^e^f-(0|V^.(0)a(P, S)a\P, 5)V'.(^)|0), (4) 

where a{P, S)[a\P, S)] is the destruction (creation) operator for the fragmented 
hadron, of given four-momentum and PL four-vector. 

The hadronic tensor (1) is not color gauge invariant. Introducing a gauge link is 
not sufficient to fulfil this condition, but EOM suggest to add suitable contributions 
of higher correlators, involving two quarks and a number of gluons, so as to construct 
a gauge invariant hadronic tensor. 

We adopt an axial gauge, obtaining for the correlator a gM/ Q expansion, where g 
is the coupling, M the rest mass of the hadron and Q the QCD "hard" energy scale, 
generally assumed equal to \J\q'^\- We examine in detail the first two terms of the ex- 
pansion. The zero order term corresponds to the QCD parton model approximation. 
As regards the second term, it concerns the T-odd functions; in particular, we discuss 
an interesting approximation, already proposed by Collins[79]. In both cases we ob- 
tain several approximate relations among "soft" functions, which survive perturbative 
QCD evolution, as a consequence of EOM. Our approach allows also to determine 
the Q-dependence of some important azimuthal asymmetries and to draw conclusions 
about the Burkhardt-Cottingham[80] and Efremov-Leader-Teryaev[65] sum rules. 

Section 2 is devoted to the gauge invariant correlator (more properly to the dis- 
tribution correlator), whose properties are deduced with the help of EOM. In partic- 
ular, we derive an expansion in powers of gM/ Q, whose terms can be interpreted as 
Feynman-Cutkosky graphs. In section 3 we give a prescription for writing a gauge 
invariant sector of the hadronic tensor which is of interest for interactions at high Q. 
In sects. 4 and 5 we study in detail the zero order term and the first order correction 
of the expansion, deducing approximate relations among functions which appear in 
the usual parameterizations of the correlator [55, 81]. Sect. 6 is dedicated to the frag- 
mentation correlator. In sect. 7 we illustrate the azimuthal asymmetries involved in 
the three different deep inelastic processes. Lastly sect. 8 is reserved to a summary 
of the main results of the paper. 
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2 Gauge Invariant Correlator 



The correlator (3) can be made gauge invariant, by inserting between the quark fields 
a link operator[76, 77, 55], in the following way: 

*,,(p;P,5) = nJ ^^e^P^P,S\i^^{0)£{x)i;,{x)\P,S). (5) 

Here 

C{x) = Pexp[igAjix)], with Ax(x) = f A„A"(^)rfz^ (6) 

is the gauge link operator, "P" denotes the path-ordered product along a given in- 
tegration contour T, and being respectively the Gell-Mann matrices and the 
gluon fields. The link operator depends on the choice of T, which has to be fixed 
so as to make a physical sense. According to previous treatments [5 5, 79, 82, 83], 
we define two different contours, T±, as sets of three pieces of straight lines, from 
the origin to Xioo = (±oo, 0, 0^), from Xioo to X2oo = (±oo, x^) and from X2od to 
X = {x^ , x^^), having adopted a frame, whose 2;-axis is taken along the hadron mo- 
mentum, with x"^ = l/^/2(t±z). We remark that the choice of the path is important 
for the so-called T-odd^ functions [4 7]: the path 2+ is suitable for DIS distribution 
functions, while T_ has to be employed in DY[82, 83]. For an antiquark the signs of 
the correlator (5) and of the four-momentum p have to be changed. 

In the following of the section we investigate some properties of the correlator. 



2.1 T-even and T-odd correlator 



We set [82] 



^E{0) = -[$+±$-], 



(7) 



where corresponds to the contour I± in eqs. (6), while $^ and $o select re- 
spectively the T-even and the T-odd "soft" functions. These two correlators contain 
respectively the link operators Ce{x) and Co{x), where 



^E{o){x) ^ {exp igKx^{x) ± exp igAi_{x) } 



(8) 



■fMore precisely, one should speak of "naive T", consisting of reversing all momenta and angular 
momenta involved in the process, without interchanging initial and final states[84, 85, 86]. 
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and ^x±{x) are defined by the second eq. (6). Eqs. (7) and (8) imply that the T-even 
functions are independent of the contour (2+ or while T-odd ones change sign 
according as to whether they are involved in DIS or in DY[79, 82]. In this sense, 
such functions are not strictly universal[79], as already stressed. It is convenient to 
consider an axial gauge, 

A- = A+ = 0, (9) 

with antisymmetric boundary conditions[55]. Here we have adopted the shorthand 
notation A'* for A"74^. In this gauge - proposed for the first time by Kugut and 
Soper[87] and named KS gauge in the following - we have 

Ax+(x) = -Ai_(x) = / dz^X^^iz), (10) 

where Xi is a shorthand notation for Xj^+oo, i = 1, 2. Therefore, in the KS gauge, 

Ce{x) = Pcos [fiiAx+(x) , Co{x) = iPsin \gKx^{x) (11) 

and the T-even (T-odd) part of the correlator consists of a series of even (odd) powers 
of each term being endowed with an even (odd) number of gluon legs. As a 
consequence, the zero order term is T-even, while the first order correction is T- 
odd. This confirms that no T-odd terms occur without interactions among partons, 
as claimed also by other authors[57, 58, 59, 79]. Gauge invariance of the correlator 
implies that these conclusions hold true in any axial gauge, such that conditions (9) 
are fulfilled. From now on we shall work in such a type of gauge [88, 89]. 

2.2 Power Expansion of the Correlator 

We consider which, as explained before, refers to DIS. As regards DY, the T-odd 
terms will change sign, as follows from the choice of the path - X_ instead of X+ - and 
from the first eq. (10) and from the second eq. (11). We rewrite C{x) as 

oo 

C{x) = Y^iwYKix). (12) 

n=0 

Here Ao(x) = 1, while for n > 1 one has, in the KS gauge, 

l'X2 fZl fZn-1 

An{x)= dz^^ dz^\.. dzi:-[AM...A,,{z2)A,^{z^)], (13) 

Jxi Jxi Jx\ 



where the Zi = (oo, z^, Zi±), i = 1,2, ...n, are points in the space-time along the hne 
through xi and X2. Substituting eq. (12) into eq. (5), we have the following expansion 
of $ in powers of g: 

00 

<^> = E(^^?rr„, (14) 

n=0 

with 

74 

(r„)., = iv| ^e^^'-(P,5|^,.(0)A„(x)^,(x)|P,5). (15) 

As already noticed, is T-even for even n and T-odd for odd n. 

Now we invoke the Politzer theorem[71], concerning EOM. This states that, if 
we consider the matrix element between two hadronic states of a given composite 
operator, constituted by quark and/or gluon fields, each such field fulfils EOM, no 
matters if the parton is off-shell and/or renormalized. We show in Appendix A that, 
owing to the Politzer [71] theorem, the term Tq fulfils the Dirac homogeneous equation, 
i. e., 

{^-m)To = 0, (16) 

where m is the quark rest mass. The corresponding Feynman-Cutkosky graph is 
represented in fig. 1. For n > 1 we have instead 



Figure 1: Feynman-Cutkosky graph for zero order term of expansion (14). 



(17) 
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Here we have set 

gm-nn ^ a. 

p-m + ie - It^ - m + ie 

X , ^ --^7'^^ (19) 

I 

ki = J2kr. (20) 

r=l 

The kr (r = 1,2, ...n) are the four-momenta of the n gluons involved in the quark- 
gluon correlator „ , defined as 



J{p-kn)x 



X (P,5|^,.(0)P'[B^„(fcO...B^,(fci)]^,(a;)|P,5), (21) 



with 



B^(A;) = A^(A;) + A^(A;), (22) 

A,(^) = / (04AM(^)e^'^ (23) 

kJk) = 6{k+) lim [ dKe-'^^AJk-,K,k^), (24) 

the reference frame being the one defined at the beginning of this section. Moreover 
the operator product P' is defined according to the following rules: 

- any A^(A;) is at the left of any A^(A;); 

- the Af,{k) are ordered as A^^(/ci)A^2(A;2)...A^;(A;,); 

- the A^(A;) are ordered as A^^{km)...A^.^{k2)A^^{ki). 

Lastly the quark-gluon correlators fulfil the following homogeneous equa- 

tion: 

(l^ - ^„ - m)$W..^Jp, /ci, k2...kn) = 0. (25) 

Each term of the expansion (14) - somewhat similar to the one obtained by 
CS[76, 77] - may be interpreted as a Feynman-Cutkosky graph. It corresponds to 
an interference term between the amplitude 

" nucleon ^ quark -I- spectator partons" (26) 



without any rescattering, and an analogous one, where n gluons are exchanged be- 
tween the active quark and the spectator partons. In particular, the interference term 
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Figure 2: Same as fig. 1 for first order correction in the coupling. 



is such that the gluons (for n > 0) are attached to the left quark leg, see figs. 2a and 
3a. An important result, deduced at the end of Appendix A, is that such a term turns 
out to correspond to any interference term between two amplitudes, such that k and 
n — k gluons respectively are exchanged between the active quark and the spectator 
partons, with < k < n. The situation is illustrated in figs. 2 and 3 for n = 1 and 2. 

It is worth noting that a radiation ordering similar to the one established here is 
found in semiinclusive processes at large a; [90] and in totally inclusive DIS at small 
x[91]. Moreover the terms (21) consist of quark-gluon-quark correlations, analogous to 
the one introduced by Efremov and Teryaev[75] and by Qiu and Sterman[50, 51, 52]. 

As a consequence of the Politzer theorem, formulae (14) to (21) hold for renormal- 
ized fields, provided we take into account the scale dependence of the coupling g, of 
the quark mass m and of the correlators ^i; ^2---^n)[92]. Moreover one has 

to observe that the four-momenta appearing in the propagators are highly off-shell: 

and {p — kr)^ are of order Q'^[77, 78], because the uncertainty principle demands 
hard interactions to occur in a very limited space-time interval, corresponding to the 
condition 

l/l > Ml (27) 
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Figure 3: Same as fig. 2 for second order correction. 



Therefore we have 



2p^p and = 0{Q), whence 
\P-\ = 0{Q) 



{21 



and it follows that the coefficients r„ are of order Q"", up to QCD corrections, 
consisting of terms of the type g'^^{lnQ)"^, with k and m integers and k > m[93]. 
For the same reason, the coupling g, which appears in expansion (14), assumes small 
values, corresponding to short distances and times. 

To summarize, we have found that the T-even and the T-odd correlators, given 
by eqs. (7), may be written as expansions in gM/Q, i. e., 

^.(p)-Ef^l^., .o(.)^f:f^V""r_, (29) 



Q 



Q 



n=0 \ / n=0 

where = F„Q"/M" depend still on Q, as told above. As explained above, $o 
changes sign when involved in DY. Stated differently, T-odd terms present an odd 
number of quark propagators, see eq. (19) for odd n: in the limit of negligible quark 
mass, quark four-momenta in DIS are spacelike, whereas in DY they are timelike[82]. 

The first two terms of expansion (14) will be studied in detail in sects. 4 and 5 
respectively. 
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3 Hadronic Tensor 



In the present section we refer indifferently to the hadronic tensor of the three pro- 
cesses introduced, which may involve one or two fragmentation correlators; in fact, 
as we shall see in sect. 6, this object requires only minor modifications with respect 
to the treatment of last section. 

If we substitute the gauge invariant correlator (5) into the hadronic tensor (1), this 
latter does not fulfil the requirement of electromagnetic gauge invariance: only the 
term of zero order in the coupling satisfies this condition. In order to get a complete 
gauge invariance at any order, we have to recall the interpretation given above of 
the correlator. For example, at first order in the coupling in SIDIS, we see that the 
"hard" scattering amplitude 97* — q'g - where we have denoted by q and q' the initial 
and final quark and by ^ a gluon - consists not only of the graph of fig. 4a, encoded 
in the first order term of the correlator, but also of the one represented in fig. 4b, 
which interferes coherently with it. This guarantees electromagnetic gauge invariance 
for the first order graph[94]. Furthermore, convoluting "hard" graphs with the "soft" 




Figure 4: Graphs for "hard" amplitudes interfering coherently, first order correction 
in the coupling. 



factors, these two amplitudes give rise, among other objects, to asymmetric Feynman- 
Cutkosky graphs (fig. 5), related to interference terms. These are observables - 
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necessarily gauge invariant - and therefore assume real values. This procedure, already 
suggested in ref. [78], can be generalized to the three kinds of hadronic tensors 
considered in the present article, at any order in g, so as to obtain sets of graphs 
corresponding to observable, and therefore gauge invariant, quantities. We show how 
to construct them at any order n, corresponding to the overall number of gluons 
exchanged between active quarks and spectator partons. The procedure consists in 
following steps, for a given n: 

- Consider the n + 1 possible combinations of gluons occurring in the hadronic 
tensor (1), say, s for hadron A and n — s for hadron B, with s = 0, 1 ... n. 

- For a given s (n — s), consider all possible correlators, according to the definition 
given in subsect. 2.2: as seen at the end of last section, these amount to s+1 (n— s+1) 
correlators equal to Tg (Tn-s)- 

- Add each such correlator all those graphs whose "hard" parts interfere coherently 
with it, as shown in fig. 5. In practice, one has to do this for the correlator whose 
gluons are attached to the "left" quark leg and to multiply by the number of gluons 
of each correlator. 




Figure 5: Feynman-Cutkosky graphs corresponding to "hard" amplitudes of fig. 4. 
Also the complex conjugate graphs, which amount to specular images of these two, 
contribute to first order corrections. 

Then we have, up to QCD corrections at each order of the expansion, 

oo 

W.M = T.wSiq), (30) 

n=0 
12 



with 

/„ 11 n 

dp-<fpj dnnY.Y.TrM^:l (31) 



r=0 5=0 



= T.i^ + ^)in-s + l)\r^:''^^frp'^^-/], (32) 



s=0 

I 



pii,r) ^ Y^STlpSl. {33) 

•m=r 

Here we have used the following shorthand notations: 

cm _ q^lr+l,|J.r+2,■■■,^^m ^n,s _ (|jn,S (0A\ 

moreover we have set SIP' — 1 ior m — r. 

Unless we integrate over some final transverse momentum [of the lepton pair in the 
case of DY, of a final hadron in SIDIS or e^e~ annihilation], the phase space of the final 
gluons emitted undergoes a restriction [93], expressed by a doubly logarithmic form 
factor; this is more and more sizable at increasing energy, resulting in the well-known 
Sudakov-like damping[76, 95]. Lastly, the complicated perturbative QCD corrections 
to the parton distribution functions may be calculated according to the algorithm 
suggested in refs. [76, 77]. 



4 Zero order term: the QCD parton model 

In this section and in the following one we shall be concerned with the hadronic 
tensor for DY process, but our results may be trivially extended to the other two 
deep inelastic processes, with a slight difference for the fragmentation function, to be 
discussed in sect. 6. 

Let us consider the hadronic tensor (31) at zero order, i. e., 

W^'^ ^Cj £p^dp-Tr [7ar^(p)7/3ro^(p')] , (35) 

where the Fo's are given by eq. (15) for n = and fulfil the homogeneous Dirac 
equation (16). The tensor (35), T-even, corresponds to the Born approximation 
considered in the introduction. As appears from eq. (35), the study of this tensor 
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amounts to analyzing the correlator Fq, which is itself T-even and gauge invariant at 
zero order in g. In Appendix B we show that 

ro(p) = ^{]^ + m) [hip) + 75^^il(p) + 75^/iiT(p)] 2p+S{p' - w?). (36) 
Here fi{p), giiip) and hi^ip) are functions of the four-momentum p of the active 



quark, which, in this case, is on shell, p = {E,p), with E = ^Jrn?^ + p^. S^^ and S]_ 
are the components of the quark PL vector, respectively parallel and perpendicular 
to the hadron momentum. Moreover we have set 

V^^p-n., (37) 

having defined the dimensionless, light-like four- vectors n± in such a way that 

n+ • n_ = 1 (38) 

and their spatial components are along (+) or opposite (-) to the hadron momentum. 
It is important to notice that, if integrated over p^, the expression obtained for the 
zero order correlator turns out to be proportional to the density matrix of a quark 
confined in a finite volume, but free of interactions with other partons[96]. Therefore 
we fix the normalization constant so as to obtain, after integration, just the density 
matrix i. e., 

N = 2r. (39) 

Lastly, it is convenient to express S'^^ and S'\_ in terms of the components of the PL 
vector of the hadron. As shown in Appendix B, one has 

Sl = \ (£) -r)^ + 0{vl), Sl^S^ + Ax£ + 0{vl). (40) 

Here 

A ^ S.^S-X^±^, (41) 

p = {\p\,Ep), P = p/|p|, rj^^pjV, (42) 

A± = -S-rJ^, p^ = (0,0,p^) (43) 

and p± is the transverse momentum of the active quark with respect to the hadron 
momentum. 
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Equation (36) has important consequences on TMD T-even functions, as we are 
going to illustrate in the two next subsections. To this end we compare that equation 
with the naive parameterization of the TMD correlator in terms of Dirac components, 
without introducing any dynamic conditions [55, 56, 81]. We give such a parameteri- 
zation in Appendix C, up to and including twist-3 terms. The twist-2, T-even sector, 
which we study in subsect. 4.1, corresponds to quark distribution functions which 
survive when interactions with gluons are turned off. As regards the twist-3 func- 
tions, we distinguish among the T-even, the T-odd and the "hybrid" ones, these lasts 
deriving contributions both from T-even and T-odd terms. 

4.1 Twist-2, T-even Correlator 

If quark-gluon interactions are neglected, the correlator is usually parameterized as[97, 
98] 

f V 1 

= Tf^-^^^^ ^ ^^^^^ ^ ^^9iT)lbi+ + 2^it75[(^±, i+] 

+ ^(A/i^^ + A^/i^^)75[?^±,^i+]}2p+7r5(p^-m2). (44) 

Here the Dirac operators considered are purely T-even, as can be checked; moreover 

^± = P±//"o, A_L = --S • r^_L (45) 

and jiQ is an undetermined energy scale, introduced for dimensional reasons, in such 
a way that all functions embodied in the parameterization of $ have the dimensions 
of a probability density. This scale [97] determines the normalization of the functions 
which depend on r]±. In particular, as is well-known, the 6 twist-2 functions, which 
appear in the parameterization (44), are interpreted as TMD probability densities: 
/i is the unpolarized quark density, qil the longitudinally polarized density in a 
longitudinally polarized (spin 1/2) hadron, qit the longitudinally polarized density 
in a transversely polarized hadron, h^^ the transversity in a longitudinally polarized 
hadron and 

KT-hiT + WAh^T (46) 
is the TMD transversity in a transversely polarized hadron. 
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Now we compare the parameterization (44) with the correlator (36). To this end 
we consider projections of both matrices over the various Dirac components, i. e., for 
a given Dirac operator F, 

$^ = -TrT^, (47) 

possibly taking into account eqs. (40). 

First of all, F = 757^ and 757"*" 7i (i = 1, 2) yield, approximately in the limit of 
m = 0, 

2 

hiL ~ -^giL, giT ~ iphiT, Kt ~ ^^12^- (48) 

These relations hold up to terms of order [gM/QY, since, as we have seen, the T- 
even Dirac components of $ derive contributions only from even powers of gM/Q. 
Moreover, the Politzer theorem implies that the relations are not modified by renor- 
malization effects, and therefore hold also taking into account QCD evolution. 

In order to determine /^o, we observe that the functions involved on both sides 
of eqs. (48) are independent of V. Therefore we must set //q = CqV, Cq being 
a dimensionless numerical constant, independent of momentum. But since these 
functions are quark densities, they should be normalized adequately, setting Cq — 1. 
Then, neglecting the quark mass, 

f^o^V^^p-n^. (49) 

This result differs from the treatments of previous authors[55, 81], who assume /iq — 
M. 

By comparing CLAS[26] and HERMES [20] results, at not too high values of Q"^ 
(1.5 to 3 GeV) the first relation (48), together with eq. (49), is verified for x < 
0.35[96], discrepancies at larger x being attributed to higher twist contributions. 

4.2 Twist-3, "Hybrid" Correlator 

Now we consider a sector of the correlator which, as already explained, has both T- 
even and T-odd contributions. In particular, here we focus on that part of "hybrid" 
correlator which comes from the so-called "kinematic" twist-3 terms. In Appendix C 
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we find, according to the usual notations[55, 81], 

+ ia;M(e + ^^75^x + i(A/ii + AxM75[?^-,^+])}2p+5(p'-m2). (50) 

Comparing the projections of the operators (36) and (50) over F = 7j (i = 1, 2) yields 
the approximate relation 

^ /i, (51) 

which corresponds to the Cahn[53, 54] effect and is approximately verified for suf- 
ficiently large and small x[45]. Also this equation, hke eq. (48), survives QCD 
evolution. As we shall se in the next section, eq. (51) holds up to terms of order 
gM/Q, since /-*" derives also (T-odd) contributions from one-gluon exchange. 

The projections of the same operators over F = 757^ (i = 1, 2) yield (after inte- 
gration over p±) 

9t{x) ~ -^hiix), (52) 

where 

grix) = j d^p±g!r{x,pl) (53) 

and 

hi{x) = j (fp±[hiT{x,pl) + \ril\hiTix,pl)], (54) 

as obtained by integrating eq. (46) over transverse momentum. In this case the 
contribution of the QCD parton model is very small: m is negligible for u- and d- 
quarks, while for s-quarks hi is presumably small, because the sea is produced mainly 
by annihilation of gluons, whose transversity is zero in a nucleon. Therefore the 
contribution of quark-gluon interactions, neglected in the approximation considered, 
becomes prevalent in this well as for F = 1 and 757+7., corresponding 

respectively to e and Hl. The effect of such interactions will be discussed in sect. 5. 

4.3 Remarks 

To conclude this section, we sketch some consequences of our theoretical results. 
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A) In the expression (46) or (54) for transversity, the second term is due to a 
relativistic effect. To illustrate this, consider a transversely polarized hadron. The 
longitudinal polarization of the quark, due in this case to the transverse momentum, 
is magnified by the boost from the quark rest frame. This additional polarization, 
along the quark momentum, has again a transverse component with respect to the 
nucleon momentum. 

B) Eq. (54), together with the last two eqs. (48), suggests a method for deter- 
mining approximately the transversity of a hadron. Indeed, qit can be conveniently 
extracted from double spin asymmetry[99, 100, 101] in SIDIS with a transversely po- 
larized target. This asymmetry is expressed as a convolution of the unknown function 
with the usual, well-known fragmentation function of the pion. Therefore the method 
appears complementary to the one usually proposed[16, 102], based on the Collins 
effect [46] in single spin SIDIS asymmetry; in this latter case one is faced with the 
convolutive product of hiT with the Collins function, which is poorly known[63, 64]. 

C) Eq. (52) estabhshes a relation between transversity and transverse spin. In- 
deed, the two quantities are related to each other. But, unlike transversity, the 
transverse spin operator is chiral even and does not commute with the free hamilto- 
nian of a quark [42]: in QCD parton model it is proportional to the quark rest mass, 
which causes chirality flip. 

D) We note that giT, h^^ h^j^ are associated with "twist-2" Dirac operators[42, 
43], and yet, in our treatment, they are multiplied by inverse powers of Q {Q~^ for the 
first two functions, Q~'^ for the third one). This would be unacceptable for common 
distribution functions, but, when transverse momentum is involved, also the orbital 
angular momentum plays a role. To illustrate this point, we recall that the quark 
distribution functions may be regarded as the absorptive parts of w-channel quark- 
hadron amplitudes [44]. For example, Qit corresponds to an amplitude of the type 
(-|- -|- I — h), denoting by |AA) a state in which the nucleon and quark helicities are, 
respectively, A and A. The amplitudes corresponding to the functions in question 
involve a change AL = 1 (for gyr and h^^) or AL = 2 (for h^j) in the orbital angular 
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momentum, therefore they are of the type 

A = A{sine)^^, (55) 

where 6 = arcsin\p±\/\p\ is the angle between the nucleon momentum and the quark 
momentum, while A is weakly energy dependent. Since |p| is of order Q and |p_l| of 
order M, our result can be understood-^-. 

5 First Order Correction 

The first order correction in g of the hadronic tensor reads [see eqs. (31) and (32)] 

= -2gC I dp-d'p^ J ^JrN^p, (56) 

■with 

N^^ = 2mp,p\ k)^%{p, k)^pT^{p') + 7ar^(p)/i^(p',P, fc)<(p', k)] (57) 

and 

KiP,P',k)=^a2 ^7^ + 7^ . / ^. 7a. (58) 

P — m + te p' — f — m + te 

Moreover the $^^^'s are given by eq. (21) for n = 1 and fulfil the homogeneous Dirac 

equation 

(^-^-m)^'j}\p,k) ^0. (59) 

In Appendix B we show that, according to the Politzer theorem, and adopting, as in 
the previous sections, an axial gauge, (p, k) is parameterized as 

(P, k) = *,(p, k)6 (p^ - ""'^P'^ ^ . (60) 

Here 

Pi = p — k, with pI = (61) 

and 

k) ~ ^(l^i + m)L[C^ + AC^75^ + A^C^ts^ + I^tC'^i^M (62) 



■l-This observation is the fruit of a stimulating discussion -with Nello Paver. 
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Moreover we have 



\l\A\S^a = ^appan^nlpl, (63) 



where Vi = p^/\/2, while Lp and a are defined in Appendix B. Lastly 

= + Wan^(C2^fVx+C3M5lX), (65) 

AC^ = ACpix^, (66) 

AtC^ = AtCpix^, (67) 

AtC; = AtC'pix^. (68) 



Here Cj (i = 1, 2, 3), AC, A^C and A^C are correlation functions, the Ci being unpo- 
larized, while the others are polarized. More precisely, AC and A^C, which are, re- 
spectively, longitudinally and transversely polarized, are related to an overall nucleon 
polarization. On the contrary, AtC is a transversely polarized correlation function 
connected to quark-gluon interaction, for example, to a spin-orbit coupling[57, 58, 59]. 

5.1 Approximate Factorization 

The second term of eq. (58) is not factorizable, in agreement with the observations of 
various authors[57, 58, 59, 103, 104], who have shown failures of universality[103, 104] 
at large tranverse momentum. However, for sufficiently large Q, and adopting an 
axial gauge, this term is negligibly small[94] in comparison with the first one, which 
instead is factorizable. In fact, the gluon corresponding to the first term has a smaller 
offshellness than the one involved in the second term. This approximation is especially 
acceptable, even for relatively small provided we limit ourselves to small transverse 
momenta [79] of the initial hadrons with respect to the direction of the momentum 
of the virtual photon in the center of mass of the DY pair. However, as already 
explained in sect. 2, also in the case when factorization is approximately satisfied, 
the T-odd distribution functions change sign from SIDIS to DY. We shall illustrate 
phenomenological implications of this change of sign in sect. 7. 
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In this approximation eq. (57) simplifies to 
with 



ri(p) 



-r 



(69) 



(70) 



m + ie' J {271^ 

Then, in an axial gauge, under the kinematic conditions above described, the tensor 
W^p can be written [see eqs. (56) and (69)] in a form similar to W^p, giving thus 
rise to an approximate [5 7] factorization of T-odd functions. Our conclusion is quite 
analogous to the one drawn by Collins [79] and presents some similarity with the 
Qiu-Sterman assumption about the quark-gluon-quark correlation functions [50]. 

Moreover, eqs. (65) to (68), together with eqs. (40), induce for Fi the following 
parameterization, at twist-3 approximation: 



ri(p) 



ie) 2 



Here 



Moreover 



nr 



Po = IPNP- 



J o 

h' 
dn 



71 



J dnCi, 

dQArC, 

d^k p~Pi 
(27r)4 2p+ 



dQ,C3, hT,o ^- J d^lArC, 
fi^ I (ifi(C2 + rAC), 



r — 



k-p^ 



V_ 



cosh(fi ± a 



sinhip 
2ip 



and 



PiP L( ) ' 
d^k^2p+d'^Pi5{pl-m'). 



(71) 
(72) 

(73) 
(74) 

(75) 

(76) 



The notations for the functions are somewhat similar to those introduced by refs. 
[55, 81]. The suffix "o" in gT,o and /it,o denotes T-odd contribution to these three 
functions, classified as "hybrid" in sect. 4, since they have T-even counterparts, see 
eq. (50). These functions are normalized coherently with their counterparts: indeed, 
if the quark is on-shell. 



[7r(p^ — + ie) 



—iS{p^ — 



(77) 
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the (— i)-factor being compensated by the i— factor in expansion (14). This constrains 
also the normahzation of the other functions included in eq. (71). Moreover, as 
already noticed in connection with correlation functions, the function h' describes a 

quark transverse polarization induced by quark-gluon interactions: this polarization, 
present also in spinless or unpolarizcd hadrons, is somewhat similar to the Boer- 
Mulders function[47], although it is twist-3 and not twist-2. 

5.2 Twist-3, T-odd correlator 

Now we compare the parameterization (71) with the naive parameterization of the 
twist-3, interaction dependent correlator. This reads, according to Appendix C, 

$^ = $^ + (78) 

where ^'jj is given by eq. (50), substituting 5{p^ — m?) by [n{p'^ — vr? -|- ie)]~^. On 
the other hand 

+ l.^A.'^i.h' + ^75[7+, 7-]p± • (79) 

Comparison between parameterization (78) and result (71), component by compo- 
nent, yields the following approximate relations: 



9 


Jo 1 


ft ~ 9li 




(80) 






~ "T,o ~ 


hT,o, 


(81) 


p' s 








(82) 




^ f-L 


^ 9t,0 ~ ^0 




(83) 



Also these equations survive QCD evolution, hke eqs. (48) and (51). 
5.3 Remarks 

A) From most of the above relations one can see that a unique function represents 
both unpolarized (/ or e) and polarized {g or h) distribution functions, independently 

22 



of the fact that the nucleon is polarized or not. This is a consequence of the spin-orbit 
couphng[57] in gluon-quark interactions. For the same reason, the parameterization 
of includes 5 independent functions and not only 3, despite the fact that it fulfils 
the homogeneous Dirac equation like Fq. 

B) Among eqs. (80) to (83), those which concern only T-odd functions hold up 
to terms of order {gM/QY. On the contrary, those which involve "hybrid" functions 
- including cq. (51) - hold up to terms of order gM/Q. Analogous approximate 
relations of this latter type have been found in ref. [105]. 

C) By integrating the correlator (71) over the transverse momentum of the quark, 
we obtain interesting results as regards twist-3 common functions. First of all, the 
fourth eq. (83) implies that e{x) derives just T-even contributions, and therefore, 
apart from the (negligible) term illustrated in the previous section, it is essentially of 
order (gM/Q)^. On the contrary, the main contributions to gr and hi are of order 
gM/Q and are T-odd, therefore they change sign according as to whether they are 
involved in DIS or DY reaction. These last predictions could be tested by confronting 
the DIS double spin asymmetry[4, 5, 10] with the DY one[106, 107]. In the case of DY 
one has to integrate over the transverse momentum of the virtual photon; moreover, 
if possible, it may be more promising to detect t'^t" pairs, whose polarization is 
perhaps less problematic to determine[108]. 

D) Lastly, it is worth noting that, unlike previous authors[47, 56, 81], we find 
that the functions related to longitudinal and transverse polarization are associated 
to the same inverse power of Q. For example, g-^ and h' describe, respectively, the 
longitudinal quark polarization in an unpolarized nucleon. Similarly, fj^ and fx are 
unpolarized quark densities in a longitudinally and transversely polarized nucleon. 
Conversely, the twist-2 T-odd functions hj;, corresponding to transverse polarization 
in an unpolarized nucleon, and the unpolarized distribution function /ir[47] in a 
tranversely polarized nucleon find no place in parameterization (71). 
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5.4 Consequences on gi and g2 

Now we examine some consequences of our results on the DIS structure functions gi {x) 
and g2{x), whose properties have been studied by various authors[65, 66, 109]. To this 
end, here, and only in this subsection, we re-introduce the flavor indices, dropped out 
in formula (1), in order to recover the usual definitions of those functions. Recalling 
that 

gUx) = g1{x) + g^2ix) {a = % d, s), (84) 

and setting 

gi{x) = Eeak"(^)+^(^)] (^ = 1,2), (85) 

a 

9tM = E^a /^VK,o(^, Pi) +^T,o(^, Pi)], (86) 
a 

eq. (52) implies 

9M + g^ix) = gT,e{x) + gT,o{x) + O(MVg'), (87) 
where [cfr. eq. (52)] 

gTM-i:el^[K{x) + hl{x)]. (88) 

Since, as discussed in subscct. 4.2, gx^e is negligibly small for a nucleon, result (87) 
is in contrast with the Burkhardt-Cottigham[80] (BC) sum rule, i. e., 

[\2(x)dx = 0. (89) 

^0 

Indeed, integrating both sides of eq. (87) between and 1, and assuming relation 
(89), implies 

/ gi{x)dxp^ I gT,o{x)dx, (90) 
JO Jo 

which is impossible, since gi{x) is a T-cven function, while g^.o is, by definition, T-odd. 
Furthermore eq. (89) implies, together with the operator product expansion [65], 

g,{x) + g2{x)= f'%i{y)+g?\ (91) 
Jx y 
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where gj- is the twist-3 contribution to g'T[65], to be identified, according to our 
results, with qt^o- Then eq. (87) would yield 

f'%i(y)^gT,e(^) + 0(MyQ% (92) 
Jx y 

which appears in contrast with data of gi{x)[l, 2, 9], enforcing arguments against the 
BC rule (See ref. [65] and articles cited therein). An experimental confirmation of the 
violation of the BC rule was found years ago in a precision measurement of 512(2^) [10]. 
Also the Efremov-Leader-Teryaev (ELT)[65] sum rule, i. e., 

[' dxx[gi{x) + 2g2{x)]^0, (93) 
Jo 

is in contrast with our result. Indeed, it gives rise, together with eq. (87), to the 
approximate relation 

/ dxxgi{x) / dx2xgTJx), (94) 
Jo Jo 

which, again, relates a T-even function to a T-odd one. 



6 Fragmentation Correlator 

The fragmentation correlator (4) can be made gauge invariant analogously to the 
distribution correlator, i. e., for a quark, 

A,j(p; P, S)^2V j ^e^^-(0|£(x)V^^.(0)a(P, S)a\P, 5)V'^(x)|0), (95) 

where C{x) is given by eq. (6). Object (95) may be treated analogously to the 
distribution correlator, according to the previous sections. Indeed, also in this case, 
for an antiquark one has to change the four-momentum from p to —p and to put a 
minus sign in front of the correlator. Moreover one has to choose the path 2+ for quark 
fragmentation from e^e~ annihilation, whereas the path X_ refers to fragmentation in 
SIDIS. The only important difference with the distribution correlator is that one has 
to take into account also the nonperturbative interactions among the final hadrons 
produced. However, as we shall see in a moment, this does not involve any change in 
the parameterization. 
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We treat only the case of pions, adopting for T-odd terms an approximation 
analogous to the one discussed in subsection 5.1, valid for small transverse momenta 
of the final hadron with respect to the fragmenting quark. Under this condition, we 

have 

A(p) = 2p+{A^^\p)5(p^-m^) + A^\p)[7r(p^-m^ + ie)]-^}, (96) 
^^^\P) = l{^ + m)D^, (97) 
A^'^ip) = 7-7+[^±^;^+l^±^^1- (98) 

Here Z)^ is the common fragmentation function of the pion, D;^ is defined as in ref. 
[55] and H' assumes the role of the Collins [46] function, describing the asymmetry of 
a pion fragmented from a transversely polarized quark. 

Final state interactions give rise to terms - for instance interference terms - which 
decrease as inverse powers of Q, independent of the nature of the interactions them- 
selves. To show this, we observe that these interactions may produce an azimuthal 
asymmetry in a pion fragmented from a transversely polarized quark[46, 110]. Anal- 
ogously to the distribution functions illustrated in remark D at subsect. 4.3, such 
an asymmetry may be regarded as the absorptive part of an amplitude of the type 
(-|-|— ), where ± denotes the helicity of the fragmenting quark. This kind of amplitude 
- a typical helicity flip one - behaves as 

{+]-) = B sine, (99) 

where S is a given function, weakly dependent on the quark momentum, due to 
perturbative QCD corrections. Analogously to cq. (55), we conclude that interference 

terms decrease at least as Q~^, just like the second term in cq (98). Our result 
agrees with the approach by Collins and Soper[76], who do not include "soft" final 
state interaction in the leading term of (almost) back-to-back fragmentation in e+e~ 
annihilation. 



26 



7 Asymmetries 



In this section we consider some important azimuthal and single spin asymmetries, 
which, as is well known, may be produced by coupling two chiral-even or two chiral- 
odd TMD distribution or fragmentation functions. More precisely, the terms of the 
hadronic tensor which give rise to asymmetries are written as convolutive products 

of two "soft" functions and depend on some azimuthal angle 0, relative to the final 
hadron (for SIDIS and e^e" annihilation), or to the final muon pair (for DY). Some 
of these asymmetries arise from the first order correction of the hadronic tensor, 
while others belong to the second order one, whose complete parameterization is not 
considered in this paper. 
A) Cahn effect 

This effect, pointed out for the first time by Cahn [53], has been exhibited by ref. 
[45] examining some SIDIS data[13, 14, 15] (see also ref. [HI])- We consider the 
asymmetry corresponding to the "product" 

Ac^f^®D, + h®Di. (100) 

This asymmetry is proportional to cos0 and decreases like Q~^. To the extent that 
/-*- and can be approximated by /i and respectively, one speaks properly of 
Cahn effect[45]: this amounts to neglecting quark-gluon interactions, see eq. (51) for 
distribution functions, an analogous equation holding for unpolarized fragmentation 
functions. This approximation is acceptable for relatively large Q and at small as 
shown by ref. [45]. However, one has to observe that both /-*- and D-^ are "hybrid" 
functions and in general their T-odd contributions cannot be neglected. 
It is worth considering also the "product" 

Ac2^f^®D^, (101) 

which generates a cos2(j) asymmetry decreasing like ^ hardly distinguishable from 
another one, arising from the "product" of two chiral-odd functions, as we shall see 
in a moment. Under the approximation just discussed, we predict a sort of "second 
order" Cahn effect. 
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B) Qiu-Sterman effect 

An important transverse single spin asymmetry is the one predicted by Qiu and 
Sterman[50, 51, 52] (QS) (see also refs. [75, 112, 113, 82]). This can be observed both 
in SIDIS and in DY, where one integrates over the transverse momentum, respectively, 
of the final hadron detected and of the final pair. This is described by the "products" 

AQs(xg'T®D^ (in SIDIS) and oc + c.c. (in DY), (102) 

the "bar" indicating the antiquark function and c.c. "charge conjugated". A similar 
effect could be observed in e^e~ annihilation, if one of the the hadrons observed 
were spinning. This asymmetry decreases like Q~^. Moreover, since g'rp is prevalently 
T-odd, while /i, fi and D^^ are T-even, the asymmetry is expected to assume an 
opposite sign in SIDIS and DY. 

C) Sivers effect 

According to our treatment, the Sivers asymmetry[48, 49] is described by the 
" product" 

Asivo^q't^D^ (in SIDIS) and oc c/^ + c.c. (in DY). (103) 

Therefore this asymmetry - detected by HERMES [20, 21] and COMPASS [23] 
experiments - is described a bit differently than in current literature [47, 64, 45]; in 
particular it results to decrease as ^. Moreover, comparing eq. (103) with eq. 
(102) shows a relation between the Sivers asymmetry and the QS asymmetry, as 
already noticed by other authors[82, 114, 115, 116, 117]. This relation is especially 
close if one adopts the approximation /"*" ~ /i, or ^ D-j^ according to the Cahn 
effect [53, 54]. In this approximation one would observe the already predicted change of 
sign[79, 62] in the asymmetry, similar to the QS effect; but if quark-gluon interactions 
- and therefore T-odd components of such functions - are not negligible, the prediction 
is not true. 

D) Collins effect and Boer-Mulders effect 

^We obtained a different result in a previous paper [61], since we had started from a parameter- 
ization which is usuaUy assumed for the correlator, but which is not in agreement with the results 
of the present paper. 
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In the framework of chiral-odd functions, an important single spin asymmetry 
is produced by combination of two transversities. In particular, single transverse 
polarization gives rise to an asymmetry described by the "product" 

AcoL oc hiT^H' (in SIDIS), or (104) 
Abm oc hiT®h' + c.c. (in DY). (105) 

This asymmetry - exhibited by HERMES [20, 21] data - is predicted to decrease hke 

We have also the cos2(p asymmetries 

AcL2 oc h'0H' (in SIDIS), or (106) 
Abm2 oc h' ^h' (in DY), or (107) 
AcL3 OC H' ® h' (in e'^e" annihilation), (108) 

which decrease like Q"^. Also the asymmetries (104) to (108) - of which Abm2 has 
been detected experimentally[27, 28, 29] - are described differently than in other 
articles[47, 63, 118]. As regards the dependence of the Boer-Mulders asymmetry, 
our prediction is supported[61] by DY data[27, 28, 29]. On the other hand, the 
dependence of the Collins and Sivers asymmetries might be tested in new planned 
experiments at higher energies[38]. 

8 Summary 

In the present paper we have studied the gauge invariant quark-quark correlator, 
which we have expanded in powers of the coupling and split into a T-even and a T- 
odd part. Working in an axial gauge, the Pohtzer theorem on EOM has allowed us to 
interpret each term of the expansion in terms of Feynman-Cutkosky graphs, involving 
higher correlators and corresponding to powers of gM/Q. We have also elaborated an 
algorithm for writing a gauge invariant sector of the hadronic tensor in deep inelastic 
processes, like SIDIS, DY and e+e" annihilation. However, in the gauge considered, 
and especially at small transverse momenta, the terms of the correlator prevail over 
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the remaining ones, which imphes approximate factorization of the various terms of 
the expansion, as we have shown exphcitly for first order correction in gM/Q. 

The zero order term and the first order correction of the expansion have been 
examined in detail. In both cases the Pohtzer theorem produces a considerable re- 
duction of independent functions with respect to the naive parameterization in terms 
of Dirac components, giving rise to approximate (up to powers of gM/Q) relations 
among "soft" functions. These relations survive QCD evolution. One such relation 
has been approximately verified against experimental data[26, 20], another one sug- 
gests a method for determining approximately transversity, while others could be 
checked in next experiments[34, 35, 37, 15]. Also an energy scale, introduced in the 
naive parameterization for dimensional reasons, has been determined in our approach, 
leading to predictions on dependence of various azimuthal asymmetries. One such 
prediction finds confirmation in unpolarized DY data[27, 28, 29]. The hierarchy of 
TMD functions in terms of inverse powers of Q is established taking into account not 
only the Dirac operators, as in the case of common functions[42, 43], but also the p±_ 
dependence, since in this case the orbital angular momentum plays a role as well as 
spin. 

Moreover a relation is found among gx-, the QS asymmetry and the Sivers asymme- 
try; in particular, both gT and the QS asymmetry are found to change sign according 
as to whether they are observed in SIDIS or in DY. Some doubts are cast, instead, on 
the predicted analogous change in the Sivers effect. We draw also some conclusions 
about the structure function g2{x), in particular against the BC and the ELT sum 
rules. 

Quark fragmentation involves "soft" interactions among final hadrons, but this 
does not imply a substantial difference with the distribution correlator. Rather, a 
caveat should be kept in mind for timelike photons, in DY and e'^e" annihilation, 
when Q approaches the energy of a vector boson resonance, like the T or the Z^. 
Since such a resonance interferes with the photon, one has to take into account its 
offshellness, quite different than Q^. This situation deserves a separate treatment. 

As a conclusion, we stress that, although other authors already proposed, years 
ago, a decomposition of the hadronic tensor in terms of Feynman-Cutkosky graphs[72, 
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73, 75, 78], our deduction, based on EOM, leads to strong constraints on the param- 
eterization of the "soft" parts of the graphs. 
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Appendix A 

We deduce a recursion formula for the terms of the expansion of the correlator. 
Our starting point is the Politzer theorem[71], which implies 

{P,S\^^{Q)C{x){ip-m)uUx)\P,S) ^Q. (A. 1) 

Here |P, S) denotes the state of a hadron (for instance, but not necessarily, a nucleon) 
with four-momentum P and PL four-vector S. ip is the quark field, of which we omit 
the color and flavor index. — dfj, — igA.^ is the covariant derivative, adopting for 
the gluon field the shorthand notation A^ for A^Aa. For the sake of simphcity, color 
and flavor indices of the quark held have been omitted. Moreover 

oo 

£(x) = ^(i^rA„(x), (A. 2) 

n=0 

where g is the strong coupling, while Ao(x) = 1. For n > 1, in an axial gauge 
A+ = A" = 0, we have, according to the notations and deflnitions of sect. 2, 

K{x) = / d^f / dz^\.. / [A^,(^i)A^,(^2)...A^„(^„)] . (A. 3) 

It is worth observing that 

d^K = A^(x2)A„_i. (A. 4) 

Substituting expansion (A. 2) into eq. (A. 1), we get 

oo 

Y^itgr {^j{0)An{x){i^ - m)uM^) - i'^j{0)An-i{x)[i4.{x)]uMx)} = 0, (A. 5) 

n=0 

with 

A_i(x) = and Ao{x) = 1. (A. 6) 
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Eq. (A. 5) is an operator equation, to be intended in weak sense: it holds when 
calculated between hadronic states. All equations of this Appendix will be of this 
type from now on. 

Looking for a perturbative solution for the correlator in powers oi g, we set each 
term of the series (A. 5) equal to zero, i. e., 

{%'<P - m)On{x) = iJ)L{x)On-i{x), (A. 7) 

where 

[On{x)h = ^j{0)An{x)Mx)- (A. 8) 

By Fourier transforming both sides of eq. (A. 7), and recalling relation (A. 4), we get 

- m)dn{p) = 11, j ^e'^^ [A^{x2)0^-i{x) + 0^_,{x)A'^{x)] , (A. 9) 

where 

^n{p) = [ ^.e^'^'Onix). (A. 10) 
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Eq. (A. 9) can be rewritten as 



- m)6r,{p) -il.ji^, [A^^{k)dn-,{p -k) + dn-i{p - k)A^ik)] , (A. 11) 
where 

A>'{k) = [ ^e^*="A''(a;), (A. 12) 

A^ik) = 5(A;+) lim [ dKe''^^ AJk' , K,k^). (A. 13) 

M— >oo J 

Eq. (A. 11) is a recursion formula for On{p), eqs. (A. 6) constituting the first steps. 
This formula implies eqs. (16) (for n — 0) and (17) (for n > 1) in the text. In 
particular, as regards eq. (17), the quantity r„ results in 

r^^N{P,S\dr,{p)\P,S), (A. 14) 

where N is a normalization constant. The operator On{p) in eq. (A. 11) corresponds 
to a graph endowed with n gluons, such that the n-th gluon leg is attached to the 
quark leg on the left side of the graph (see figs. 2a and 3a). 
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Taking into account the hermitian character of A''(A;)k and the relation = 
7o^n(p)7o, eq. (A. 11) imphes 

On{p){i> -m)^-ij - k)k^{k) + k^\k)dr,-,{p - k)]^,. (A. 15) 

In this case On{p) corresponds again to a graph with n gluons, but such that the n-th 
gluon is attached to the quark leg on the right side of the graph. This last result 
implies that r„ represents any graph with n gluons, each gluon leg being attached to 
the left or right quark leg. 



Appendix B 

Here we deduce the parameterizations of the quark-quark correlator at zero order 
and of the quark-gluon-quark correlation, arising from first order correction. 
B.l The Zero Order Quctrk-Qucirk Correlator 

The matrix Tq{p), defind by 

fulfils the homogeneous Dirac equation 

(|^-m)ro(p) = 0, (B. 2) 

where m is the rest mass of the quark. As shown in Appendix A, this is a consequence 
of the Politzer theorem, which implies, at zero order in the coupling, 

(^-m)V(x) = 0. (B. 3) 

Therefore, in the approximation considered, the quark can be treated as if it were on 
shell (see also ref. [119])- Then, initially, we consider the Fourier expansion of the 
unrenormalized field of an on-shell quark, i. e., 

(fp 1 

where s — ±1/2 is the spin component of the quark along a given direction in the 
quark rest frame, u its four-spinor, c the destruction operator for the flavor considered 
and 



^^^^ = / (2^;i?'""^^^^^^'^^^^' 



d'p = d'p 6 Ip- - "^^^f"- ] , V = P+/V2. (B. 5) 
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As regards the normalization of Ug and Cg, we assume 



UgUs = 2m, {P, S\cl{p')cs{p)\P, S) = {2nfS^f>' - p)g,(p), (B. 6) 



and Qsip) is the probabihty density to find a quark with spin component s and four- 
momentum p = (p~,p), with p~ — (m^ + p^)/2p+. For an antiquark the definition 
is analogous, except that, in the Fourier expansion (B. 4), we have to substitute 
the destruction operators with the creation operators dl and p with —p in the 
exponential. 

Choosing the quantization axis along the hadron momentum P in the laboratory 
frame, and substituting eq. (B. 4) into eq. (B. 1), we get 



where 



(B. 7) 




X 





But owing to the second eq. (B. 6) we have 




(B. 9) 



where 



rS(p) 



N 





2V 



s 




s 




Firstly we elaborate Fg. We have 




Here S^^ is a four-vector such that, in the quark rest frame, = (0, A/|A|P), A = S-P, 
P = P/|P| and S is the unit spin vector of the hadron in its rest frame. Therefore 
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where 

h{p) = J2{P,S\cl{p)c,(p)\P,S) (B. 14) 

s 

is the unpolarized transverse momentum distribution of the quark, while 

A'q{p) =J2MP,S\cl{p)c,{p)\P,S). (B. 15) 

s 

According to transformation properties of one-particle states under rotations, one has 

B B 
|P, S) = cos-\P, +) + i|P, -)sm-, (B. 16) 

where ± denotes the (positive or negative) helicity of the hadron and B the angle 
between P and S. Substituting eq. (B. 16) into eq. (B. 15), and taking into account 
parity conservation, we get 

A'q{p)^cosBgiL{p). (B. 17) 

Here 

9il{p) = E MP, +\cl{p)cs{p)\P, +) = - E MP, -\clip)cs{p)\P, -). (B. 18) 

s s 

is the longitudinally polarized TMD distribution of the quark, the last equality fol- 
lowing from parity conservation. 

Now we consider Fq. Eq. (B. 16) yields, for B = n/2, 

\ni))-^{\+)±i\-)), (B.19) 

where |±) and | T (i)) denote quark states with spin components, respectively, along 
P and along 

= S - AP. (B. 20) 

Substituting eqs. (B. 16) and (B. 19) into eq. (B. 11), and taking into account again 
parity conservation, we get 

niP) = ^lstnBhMp)i\ T)(T I - I I), (B. 21) 

where 

hMp) = (P,-|cUp)c_(p)|P,+) = {P,+\cl{p)c+{p)\P,-) (B. 22) 
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is the TMD transversity of the quark. Returning to the Dirac notation, we have 

I TXT I = \{Hm) (l + 75^J , I i)a I = \{Hm) (l -75^J , (B. 23) 

where is such that 5*^ = (0, n) in the quark rest frame and 



n = tI^- (B. 24) 



Then eq. (B. 21) goes over into 



ro(p; P. s) = ^\sine/^Tqmi> + Hts^x- (b. 25) 

Substituting eqs. (B. 13), (B. 17) and (B. 25) into eq. (B. 9) yields 



ro -^\ii> + rn) [f, + g,,lsq + /i^^Ts^J S (p- - ^^^^^ ) , (B. 26) 



2P2 



2p-* 



having set 5J = S^^cosO and = S^sinO. Eq. (B. 26) is a solution to eq. (B. 2), 
which is a consequence of the Politzer theorem at zero order in g. Since this equation 
survives renormalization - which generally implies only a weak (5-dependence[70, 93] 
- the structure of Fq is not changed by QCD evolution. 

Lastly we deduce the expressions of the four-vectors S'^^ and in the frame where 
the quark momentum is p. In the quark rest frame we have 

5f = (0,AP), 5l^(0,Sx). (B. 27) 

In view of the Lorentz boost, it is convenient to further decompose AP and S_l into 
components parallel and perpendicular to the quark momentum. We have 



i- — 7 

IpI 

S± = Axp + Si, i:±= \S±\cosp{cospn- sinpk), (B. 29) 



AP = Acosap + Sy, S|| = -cosaj^ + sin^aP, (B. 28) 



where 



P = T^, k = nx-?^, (B.30) 

|P| |P X n| 



a 



arccos{P ■ p) and (3 — arcsin{n ■ p). (B. 31) 
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The boost which transforms the four-momentum of the quark from (m, 0) to {E, p), 
with E — y/m? + p^, changes only the components along p of AP and of S_l. In 
particular, the boost transforms the four-vector (0, p) to p/m, with p = (|p|,£^p). 
Therefore, since a and P are 0(|pj_|/|p|) and |p|/P = 0(1), eqs. (B. 27) go over into 

^\\^^{^- + )' ^1 = + + oivl), (B. 32) 

where rjj_ — p±/V and = —S ■ rjj^. 

B.2 The Quark-Gluon-Quark Correlator 

Now we deduce a parameterization for the quark-gluon-quark correlator, defined 

by 

[$W(P,A;)].. = 7V| ^e^M^(P,5|V^,(0)[A,(fc) + A,(A;)]V'.(x)|P,5). (B. 33) 

As shown in Appendix A, the Politzer theorem implies, at order 1 in the coupling, 

(^^-^-m)$«(p,^)=0, (B. 34) 

which holds also after renormalization. Therefore our line of reasoning is the same 
as for To, that is, we start from unrenormalized fields and we take on-shell quarks, 
whose field satisfies expansion (B. 4). Substituting this expansion into eq. (B. 33), 
we get 

^i\P:k) = <if,{p,k)s(^p^- '^'^J'^ y (B.35) 
'^,ip,k) = iv/|^^-A_^^^_^,^(p',^)„^(p^)^^,(p'). (B.36) 

Here d^p' and V are defined analogously to eqs. (B. 5), 

Pi^p-K V,^pt/V2 (B. 37) 

and 

As,s'Ap', k) = (P, S\c\{p')[k,{k) + k^{k)]cs>{p,)\P, S). (B. 38) 

Moreover the matrix element (B. 38) fulfils a relation of the type 

As,s'Ap'^ k) = {27rf Cs^s'Ap'^ k)5%p' - Pi - k), (B. 39) 
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where Cs,s',n{p', k) is a quark-gluon correlator and p', pi and k are defined by eq. (B. 
7). Then eq. (B. 36) yields 

N 



and 



Po = {Po,P), 

We rewrite eq. (B. 40) as 



Po 



p5_ + 



2p^ 



N 



where 



*M ^ Y.^^,sAPu^)'^s{Pi)us{po), 

s 

*i = ^Cs-sAPl^k)'^s(Pl)u^s(Po)- 



(B. 40) 



(B. 41) 



(B. 42) 



(B. 43) 
(B. 44) 



Taking into account the appropriate Lorentz transformations for the spinors involved, 
we have 



Us{pi)Usipo) 
Us{pi)U-s{po) 



2 
1 



^i + m)C/(pi,po)(l + 2s75^ll), 



:B. 45) 



ii + m)U{pi,po)'y5{cosx^o± + sinx^±). (B. 46) 



Here 



Uipi,Po) = exp 



;(0iPi - 0oPo) • a 



(B. 47) 



01 



In 



El + |pi| 



m 



Pi = (Pi±, -Pi)), 



Pi 



Pi 

IpiI 



(B. 48) 



Ei^^pj + m\ (B. 49) 



analogous definitions holding for 4>q and po. Moreover Sq^^ and S'qj^ refer to the PL 
vector of a quark with four- momentum po, directly connected with nucleon polar- 
ization; they can be related to the nucleon longitudinal and transverse PL vectors, 
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using the formulae elaborated at the end of sect. Bl. S± refers to the spin caused by 
spin-orbit coupling, 

\/\pl±\S±a = eai3jpn'lnlp^_^. (B. 50) 

Lastly, X is a real parameter. 

We assume 9o,9i « 1, where 9o and 9i are, respectively, the angle between po 
and P and the one between pi and P. Then 

U{pi,Po) ~ coship + -^70(730 + lir\)sinhip, (B. 51) 

with 

V = ^V('^o-0i)' + ^Vo0i, 9 = 9^-9^, (B. 52) 

a = 01 - 00 - - 0o6'o)> r± = ■nn'Pix - -rn"Po±, (B. 53) 

2 |Pi| IPoI 

Then results in 

-^M^ k) ~ ^(J^i + m)C[C, + AC^75^ + AtC^75^^ + AtC;75^^], (B. 54) 

where 

N 1 
^ " 2^p^ [^^^^'^ ^ 2^^°^^^" ^ 7i'^l)s«^M (B. 55) 

and 

C^iPuk) = $^C,,,,^(pi,A;), (B. 56) 

AC^{pi,k) = Y.'^sCs,sAPi,k), (B. 57) 

s 

ArC^iPuk) = $^cosxC.,-.,;.(pi,A;), (B. 58) 

AtC;(pi,A;) = ^sznxC,,_,,^(pi,A;) (B. 59) 

s 

are correlation functions. We parameterize these functions with the available trans- 
verse four-vectors, since we operate within an axial gauge: 

= Cim^ + e^,,,n^(C25fp?x+C3M5>^), (B. 60) 

AC^ = ACpix^, (B. 61) 

AtC^ = AtCpix^, (B. 62) 

AtC; = AtC'pix,.. (B. 63) 
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Here Ci, C2, C3, AC, AtC and AtC are "soft" functions of p and k. The parame- 
terization of is obtained by inserting eqs. (B. 54) and (B. 60) to (B. 63) into 
eq. (B. 35). Again, as in the case of Fq, the Pohtzer theorem, of which eq. (B. 34) 
is a consequence, imphes that renormahzation effects preserve the structure of that 
parameterization . 



Appendix C 

Here we consider the parameterization of the correlator in terms of the Dirac 
components, up to and including twist-3 terms. This parameterization is similar to 
the usual ones[47, 81], also as regards notations, except for an energy scale fiQ, which 
we leave undetermined here, and for the twist-2, T-odd sector, which we omit because 
it has no place in our procedure. The scale fiQ, usually set equal to the rest mass of 
the hadron[47, 81], is determined in the sects. 4 and 5. 

The parameterization reads 



where 



(C. 1) 



f V 1 
* = ^^-^1^+ + ^^^^^ + >^±9iTh5V^+ + 2^ir75[)^±, 

+ li^hi^ + ^±hiTM^±,i+]}, (C. 2) 

+ ^xM (^e + (7^75)&± + ^{XhL + \±hTh5[i-,i+]^}, (C. 3) 

+ l,^±Ah' + ^75[7+, 7- W • S^h'^}. (C. 4) 

Here denotes the "hybrid" term, both interaction free (^^j-, T-evcn) and inter- 
action dependent (^^, T-odd): the two terms have the same parameterization, but 
behave quite differently. For the "soft" functions we have adopted notations similar 
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to those employed in ref. [81] . Note, however, that in the expression of the func- 
tions fx, 1 g't-i h' and /i'-*- do not appear in the parameterization proposed by ref. 
[81]; on the contrary, we have not taken into account the functions h and /j^, defined 
in that reference. 
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